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Optical bandedge of diluted magnetic semiconductors: 
difference between II- VI and III-V-based DMSs 
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Applying the dynamical coherent potential approximation to a simple model, we have theoretically 
studied the behavior of the optical bandedge in diluted magnetic semiconductors (DMSs). For 
Ai I _ a: Mn :E _B VI -type DMS, the present study reveals that the linear relationship between exchange- 
spitting AE ex and the averaged magnetization \x{S z )\ widely holds for different values of x. The 
ratio, AEex / x(S z ) , however, depends not only the exchange strength but also the band offset. 
Furthermore, the present study reveals that in the low dilution of Gai-^Mn^As the optical bandedge 
exists not at the bandedge of the impurity band but near the bottom of host band. The optical 
bandedge behaves as if the exchange interaction is ferromagnetic although the antiferromagnetic 
exchange interaction actually operates at Mn site. We conclude that the spin-dependent shift of 
the carrier states between the impurity band and host band accompanying with the change of 
magnetization causes the apparently ferromagnetic behavior of the optical bandedge which was 
reported in the magnetoreflection measurement of Gai-^Mn^As. 

PACS numbers: 75.50.Pp, 71.23.-k, 71.70.Gm 
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I. INTRODUCTION 



In order to study the effect of the sp-d exchange inter- 
action between a carrier (an s conduction electron or p 
hole) and localized magnetic moments (d spins) together 
with magnetic and chemical disorder in DMSs, we have 
previouly introduced a simple model for Ai^ x Mn x B-typc 
DMSs p| . In this model, the local potentials of nonmag- 
netic (A) ions in a semiconducting compound (AB) are 
substituted randomly, with mole fraction (x), by the local 
potentials that include the exchange interaction between 
a carrier and the localized spin moment on a Mn (de- 
noted by M) ion. Thus, the potential to which a carrier 
is subjected at a site differs depending on whether the 
site is occupied by an A ion or M ion. The Hamiltonian 
H is given by 



H — ^ K £mn a lnfi a nfj, ~t~ ^ / ■ 



(1.1) 



where u n is either u A (on the A site) or u^ 1 (on the M 
site) depending on the ion species occupying the n site: 



E M a ln a n» ~ 1 E ■ s n«™ -C 1 - 3 ) 



,M 



Here, a^ nfl and a nfl are, respectively, the creation and 
annihilation operators for a carrier with spin fi on the 
n site. The transfer-matrix element between m and n, 



£ m n, is assumed to be independent of the types of consti- 
tutional atoms which occupy the m and n sites. In II- VI- 
based DMSs of the Af_ x Mn x B vl type, E A (E M ) repre- 
sents a nonmagnetic local potential on the A 2+ (Mn 2+ ) 
sites. In III-V-based DMSs such as Gai-^Mn^As, the 
spin-independent potential Em(< 0) can be regarded as 
a screened Coulomb attractive potential between a car- 
rier (hole) and the Mn 2+ ion (acceptor center). The ex- 
change interaction between the carrier and localized spin 
S„ of the Mn site n is expressed by —la^a^ ■ S n a nv , 
where <r ^„ represents the element of the Pauli spin matri- 
ces. Throughout this article, we disregard the electron- 
electron, hole-hole, and/or electron-hole interactions. 

The virtual crystal approximation (VCA) is widely em- 
ployed to describe the extended states in II-VI-based 
DMSs 2]. The VCA is a first-order perturbation the- 
ory with respect to the sp-d exchange interaction and/or 
the band offset energy (Em — Ea)- In standard VCA, 
first molecular field approximation (MFA) is applied, re- 
placing S„ by the thermal average (S„) taken over all 
Mn site. Second, the local potential u n is replaced by 
configuration averaged one, 

u v n CA = (l-x)u A +x{u™), (1.4) 

= E^ 1 ~ X ) Ea + X ( Em ~ Ia w( S *))] a nii. a W 

(1.5) 

for every site, where S z is the z-component of S„. Note 
(S x ) = (S y ) = because the magnetization is assumed 
to be along the z-axis. A major advantage of the VCA 
is to employ the periodic Hamiltonian 



H 



VCA 



E 

m,n,/i 



E* 



VCA 



(1.6) 
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instead of Eq. (|1.1|) . In the VCA picture, therefore, the 
carrier "sees" an effective potential u VCA = (1 — x)Ea + 



2 



x(Em T I(Sz)) at all site, where — (+) is for up- (down- 
) spin of the carrier. No spin-flip process of the carrier 
is considered. The VCA leads to the exchange energy 
splitting 

AE v ex CA = 2x1 (S z ). (1.7) 

Thus, the VCA apparently explains the energy splitting 
between <r + and o~ transition of A exciton in DMS 

AE = N o {a-0)x(S z ) (1.8) 

if we assume 21 to be the exchange constant Noa for 
conduction electrons and N f3 for valence electrons, re- 
spectively. Furthermore, the VCA explains the linear 
interpolation expression of the energy gap E g 

E g (x) = E g (0) + xV cS , (1.9) 

where V eS = E g (l)-E g (0). Equation (ESJ is experimen- 
tally observed for a given composition and temperature 
in most of II-VI-based DMS when the atoms of the group 
II element are replaced by Mn 

There exist, however, some experimental facts that in- 
dicate that the application of the VCA is limited. Among 
the them, we have already tackled some problems: the 
anomalous behavior of E g (bowing effect) in wide-gap 
DMSs |l|, the enhancement of Nq/3 with x — > ob- 
served in Cdi^^Mn^S U, and asymmetric splitting of 
Zeeman energy components In the previous works, 
applying the coherent potential approximation (CPA) 
to the present model, we studied the carrier states of 
paramagnetic ((S z ) = 0) and completely ferromagnetic 
((S z ) — S) cases. In the present work, extending the 
previous approach for the case of a finite (S z ), we study 
the nature and property of carrier states in DMSs. 

According to the VCA the exchange splitting energy 
is proportional to x(S z ) and the coefficient is equal to 
N (a - 13) (see Eq. JTJJ)). In Zni^Mn^Te and 
Cdi-^Mn^Te 7], however, it is reported that the spin 
splitting energy is proportional to x(S z ) but the coeffi- 
cient decreases with the increase in x. This may be ex- 
plained by taking the higher order effect of the exchange 
interaction into account, although imperfect treatment 
based on the second order perturbation was already done 

Another strong motivation of the present work is the 
elucidation of the sign and amplitude of the p-d ex- 
change interaction in Gai-^Mn^As. In the early stage 
of research, the /ferromagnetic (FM) coupling (No(3 > 0) 
was reported on the basis of the polarized magnetoreflec- 
tion measurement [9(. However, the exchange interaction 
between p holes and d spins is experimentall y p roved 
later to be crata/erromagnetic (AFM) [13, IH O- On 
the mechanism of the so-called carrier-induced ferromag- 
netism in Gai-^Mn^As, which has attracted much at- 
tention in recent years |l3> we have already proposed a 
theory based on the present model • According to the 
theory, the occurrence of ferromagnetism closely relates 
with the magnetic impurity band which is formed due 



to the incorporation of Mn to GaAs. When an impurity 
band exists, the optically observed bandedge may be dif- 
ferent from the bandedge of the impurity band. Thus, it 
is highly desirable to clarify what the optical measure- 
ment has detected. Throughout the present work, we in- 
vestigate the difference in the optical bandedge between 
II-VI-based and III-V-based DMSs. 

This paper is organized as follows. In Sec. II we 
briefly formulated the dynamical CPA on the bases of 
the multiple-scattering theory. The results and discus- 
sion for the behavior of the optically observed bandedge 
and the carrier states in II-VI-based DMS (or the case 
with no magnetic impurity level) is given in Sec. III. Sec- 
tion IV is devoted for Gai-xMn^As. Section V contains 
the conclusion. 



II. DYNAMICAL COHERENT POTENTIAL 
APPROACH (CPA) 

A. Dynamical CPA condition 

We shall confine our discussion to the so-called one- 
particle picture. Hereafter we assume the carrier is p 
hole, although the result does not depend on the charac- 
ter of the carrier. A carrier moving in a DMS is subjected 
to disordered potentials which arise not only from sub- 
stitutional disorder but also from thermal fluctuation of 
d spins through the p-d exchange interaction. Further- 
more, when magnetization arises, the effective potential 
for the carrier differs according to the orientation of the 
carrier spin. In the dynamical CPA [l^, the disordered 
potential is considered in terms of the spin-dependent ef- 
fective medium where a carrier is subject to a coherent 
potential, Ej or Ej, according to the orientation of its 
spin. The coherent potential £j (EjJ is an energy (in- 
dependent complex potential. Then, a carrier moving in 
this effective medium is described by the unperturbed 
Hamiltonian K: 

K = $>* + E /*)°k a *M • (2.1) 

kfi 

Thus, the perturbation term V(= H — K) is written as 
a sum over each lattice site: 

V = X>« . ( 2 - 2 ) 

n 

where v n is either v A or v^f , depending on the ion species 
occupying the n site: 

v n = ^2( E A - S^a^a^ , (2.3) 

(2.4) 
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Next, using the reference Green's function P given by 

1 



P{u) = 



u> — K 



we define the matrix t A which represents the multiple 
scattering of carriers due to the A ion potential embedded 
in the effective medium by 



f: 



v*[l-Pv. 



A-i-l 



(2.6) 



and the matrix t M which represents the multiple scatter- 
ing of carriers due to the M ion potential embedded in 
the effective medium by 



t 



M 



vt'W-Pv. 



Mi-1 



(2.7) 



Note that K, and thus P, includes no localized spin op- 
erator, and that t A (t$f) represents the complete scatter- 
ing associated with the isolated potential v A (v^f) in the 
effective medium. According to the multiple-scattering 
theory 0,0], the total scattering operator T, which is 
related to G = 1 / (u> — H) as 



G = P + PTP , 
is expressed as the multiple-scattering series, 



(2. 



n n m (^n) 



(2.9) 



n m (#n) 

Within the single-site approximation, the condition 

(in)av = for any site n (2-10) 

leads to (T) av = and thus (G) av = P. Here, we ex- 
press the average of t n over the disorder in the system 
as (t n ) av . Since the present system includes both sub- 
stitutional disorder and the thermal fluctuation of the 
localized spin an M site, the average of the t matrix is 
written as 



(l-x)t A + x(t™) 



(2.11) 



Here, (1 — x) and x are the mole fractions of A and 
M atoms, respectively; (t M ) means the thermal average 
of t over fluctuating localized spin. In the dynami- 
cal CPA, the coherent potential S M is decided such that 
the effective scattering of a carrier at the chosen site 
embedded in the effective medium is zero on average. 
Note that the thermal average of off-diagonal t-matrix 
elements (tM) = (tH) = because the magnetization is 
assumed to be along the z-axis. Therefore, the dynamical 
CPA condition is given by 



(1 
(1 



(2.12a) 
(2.12b) 



For simplicity, the t matrix elements in the site rep- 
resentation < n[i\t\nv > (n is a site index, /i, v =| or 



(2 5) I) are wr hten as t M „. The explicit expressions for t A v 



and t^f v are given in Appendix A. It is worth noting 
that in the expression of tM (tM), the spin- flip processes 
are properly taken into account. As a result, a single t- 
matrix element t*^ depends on b oth St and Sj,. There- 
fore, we solve Eqs. <|2.12a[l and (|2.12b|) simultaneously. 
Note that the diagonal matrix element tM, involves an 
operator S z which takes the quantum values of 2S + 1; 
S z = —S,—S+l,---S. Thus, the thermal average over 
the fluctuating localized spin is taken as 



S 

(t M ) = V t M 

V MA" MM 

S z = -S 



hS z 



I exp 



hS z 



a / S r =-S v ' 

(2.13) 

where h denotes the effective field felt by the localized 
spins. Since there is a one-to-one correspondence be- 
tween (S z ) and the parameter A = h/ksT 7 we can de- 
scribe the carrier states in terms of (S z ) instead of A. 
In this work, we treat the localized spins classically for 
simplicity; the actual calculations were performed for 
S = 400. 



B. DOS and local DOS 

Throughout this work, we assume the model density 
of states of the semicircular form with a half-bandwidth 
A, 



(2.14) 



as an undisturbed density of states. Then, the density of 
states with /i spin, D^(lu), is calculated by 



DM 



1. 



Im / de 



Pie) 



7r J_ A cj-e-E M (w) 



(2.15) 



for the S M determined by the CPA. In all of the present 
numerical results, we have numerically verified 



Di(uj)duj = 1 



(2.16) 



The species-resolved DOS shall help us to understand 
the nature of the carrier states. We calculate the A- 
and M-site components of the DOS, (1 — x)D A (u>) and 
xDff(uj) (fj, =t or |), where D A (w) [Df(u)\ represents 
the local DOS associated with the A (M) site (see Ap- 
pendix B). Note that 

DM = (l-x)D$(u)+xDff(u) . (2.17) 

Since D A (uS) and D m {lo) are normalized, the total num- 
ber of A-site states and that of M-site states are 1 — x 
and x, respectively. 
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C. Optical absorption spectrum 

A l i_ x Mn x B vl -type DMSs are direct-gap semiconduc- 
tors, with the band extrema occurring at the T point 
@- Upon calculating the optical absorption spectrum, 
we assume that the transition dipole moments of the A 
and M ions are same. Under this assumption, the optical 
absorption spectrum is given by the k = components of 
the DOS. Since the explicit k dependence of £& is not em- 
ployed in the present framework, we assume that k = 
corresponds to the minimum point of the model band. 
Therefore, taking £n_= —A, we define the optical absorp- 
tion spectrum by: [Tg 



A^uj) = -— Im 



(2.18) 



GaAs is a direct-gap semiconductor, whose band ex- 
trema exist at the T point |19L l20| . In Gax_a;Mn x As, the 
bottom of the conduction band is still at T point [2lj . 
even though an impurity band forms above the top of the 
valence band. Therefore, we apply the optical absorption 
spectrum defined by Eq. (|2.18|) . (see later discussion in 
Sec. VI). 



Further, we calculate the spin-coupling strength Qiui) 
defined by (see Appendix B) 



Q(w) = - 



(S(u-H)a-S)/S 



(S(co - H)) 



(2.20) 



at A'i-site 



Note that Q(u>) corresponds to —(cos 9), where 9 is the 
angle between the carrier spin and localized spin at the 
M site. Therefore, Q(u>) represents the strength of the 
spin coupling at M-site. 



III. RESULTS AND DISCUSSION FOR 
II-VI-BASED DMS 



Approximate expression for the bandedge 
energy shift 



D. Optical carrier spin polarization V(uj) and 
spin-coupling strength Q{uS) 

In order to investigate the manner of coupling between 
the carrier spin and the localized spin in DMS, we define 
the optical carrier spin polarization, V(u>), by 



V{d) = 



D^ + D^lo) 



Ej, 
A 



(2.19) 



Here we briefly summarize the approximate treatment 
to consider the behavior of the bandedge energy, although 
we can numerically solve the equation for the bandedge 
energy shift. Hereafter, we set Ea = as the origin of 
the energy. For the energy of the bottom of the band, uib, 
we assume that = — A + E(o;&) in the dynamical CPA 
condition Eq. (|2.12|) . Then, when (S z ) = S, we obtain 
the approximate expression for the bandedge energy shift 
E 6 [=E(w 6 )]: ! 




(3.1) 



with Eb = Em — IS (Eb = Em + IS) for a carrier with 
t (J.) spin. On the other hand, when (S z ) — 0, we obtain 



the cubic equation for "Et' Q 



Eg - Ht(2E M + A) + E 6 < E M - IS 



A 



E M + IS ■ 



-E M x 



~x j {E M - IS) (E M + IS) + jE M 



= 



(3.2) 



r 



As long as an impurity level (band) does not appears, Eqs. (|3.1|) and (|3.2|) give good approximate values, as 
the energy shift of the bandedges calculated by using 
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shown in the present work. Hence, on the basis of the 
approximate expression we investigate the relationship 
between the exchange integral N f3 obtained by optically 
measurement and the parameters Em and IS/ A used in 
the present model. We define No/3 by using the exchange 
energy splitting at (S z ) — S, as: 



N (3 



S b (+)-E b (-) 
xS 



(3.3) 



where and ) are the solutions of Eq. I|3.1|l for 

Eb = Em + IS and Eb = Em — IS, respectively. There- 
fore Nq(3 is a function of x, and has following limiting 
values: 



N Q /3 



21 



(i+2^r-( 2 



when x — > 0, 



(3.4a) 



N (3 - 21 



when x = 1 



(3.4b) 



B. Results for IS/ A = -0.4 and Em /A = 

To our knowledge, no impurity (acceptor) level has 
experimentally observed in A^L^Mn^B^-type DMSs. 
This implies that the fitting parameters in the present 
model should be taken as \IS + E M \ < 0.5 A for II- VI- 
based DMSs. In Figs. ~ 01 we show the results for 
IS = — 0.4A and Em = 0.0 as a sampling case of II- VI- 
based DMSs. We should note here that these parameters 
may be fit to describe a DMS like Cdi-^Mn^S which 
has the rather strong exchange interaction among II- VI- 
based DMSs (see Ref. 1). 

Figure H^a) shows how the carrier band is spin- 
polarized with the development of magnetization. Note 
that even when (S z ) = 0, the band is not same as the 
model band. Owing to the disorder of random distri- 
bution of M ion and fluctuation of localized spins, the 
band has already broaden and the bottom of the band 
has shifted to lower side from Ub = — A. With increase in 
(S z ), the bottom of down-spin band shifts to lower-energy 
side while accompanying the energy shift of bottom of the 
up-spin band. The both bandedges agree with each other 
except the case of (S z ) = S, although the down-spin 
band is strongly suppressed in the band tail. This is be- 
cause the spin-flip of up-spin carrier occurs at the energies 
wherein Dy (uj) takes a finite value. In magnetoabsorption 
and/or magnetoreflectivity spectra, on the contrary, the 
spin splitting band is observed. Thus, the present result 
for the bandedge shift is very different from the behavior 
expected from magneto-optical measurements. In Fig. 
(b), we show the optical absorption spectrum A(u>), that 
is related to the dipole transition at T point. A peak 
is found in the up- and down-optical absorption spec- 
trum, respectively. Hence, we regard the energy, w p (up) 



II-VI-based DMS x = 0.05 

IS / A = -0.4 E M / A = 0.0 
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FIG. 1: Results with x = 0.05 as a function of to/ A for 
various values of magnetization: (a) density of states D(u) 
(b) optical absorption band A(lu) in arbitrary unit (arb. units). 
Solid line is for down-spin carrier, and dotted line is for up- 
spin carrier. The approximate values of the bandedge energy, 
uJb/A, calculated using Eqs. lEnj and EH), are dotted on the 
line of {S z ) = S and (S z ) = 0, respectively. Note that the 
energy of the bottom of the model band is u — —A. 



[ojp(down)], at which the up (down)- spectrum takes a 
peak, as the up (down)- bandedge energy experimentally 
observed in optical measurement. In Fig. [2 we display 
the optical bandedge energies, w p (up) and o; p (down), as 
a function of the (S z ). The result shows that w p (up) and 
Wp(down) are roughly linear in (S z ). The behavior of the 
optical bandedge energy reproduces well an asymmetrical 
splitting of Zeeman energy component; when magnetiza- 
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FIG. 2: Optical bandedge energies lu p /A as a function of 
(S z )/S. Solid line is for down-spin carrier, and dotted line is 
for up-spin carrier. Error bar represents a half-peak width. 
The approximate values of W(,/A, calculated by Eqs. 13. II 
and 13.21 . are dotted on the line of {S z } = S and (S z ) = 0, 
respectively. 



FIG. 3: Exchange splitting [cj p (up) — (i> p (down)]/A as a func- 
tion of x{S z )/S for various values of x. The straight lines are 
adjusted to the best fit the each x data. The straight lines 
of x = and x — 1 exhibit the limiting cases given by Eqs. 
COi) and ll3~4bTl . 



tion arises the energy-shift pattern a + and a~ transition 
term of the A exciton is asymmetric relative to the posi- 
tion of (S z ) = 0. The asymmetric splitting of Zeeman en- 
ergy component has observed not only Cdi-^Mn^S 
but also in Cdi-^Co^Te j2^|. The half-peak width in- 
dicated by error bars in Fig. [21 shows that the peak of 
the optical spectrum with up-spin is broader than that 
of down-spin, which may explain the reason why the <j~ 
peak is broader than a + peak in the magnetoabsorption 
spectra 

The VCA presents the picture that a carrier in a DMS 
moves freely in an uniform medium of effective poten- 
tial u VCA — ^fxI(S z }; the depends on the orientation of 
the carrier spin. The present study reveals the feature 
of the carrier in II-VI-based DMS as below. A carrier 
does not stay with the same probability on each site, but 
tends to reside longer at Mn site due to the exchange 
interaction; Although x — 0.05 is assumed, the ratio of 
the Mn-site component of the DOS to the total DOS 
R(u>) = 0.2 ~ 0.3 at the band tail, as shown in Fig. 
??(c). The carrier state at A-site shows similar (S z ) de- 
pendence to that at Mn site (compare Figs. ??(a) and 
(b)). The result of Q(w) w 1 (see Fig. ??(b)) exhibits 
the strong antiparallel spin coupling (AP-coupling) be- 
tween carrier's spin and the localized spin realizes at Mn 
site, and the high values of the optical spin polarization 
V(u>) (see Fig. ??(a)) suggests that the carrier itinerates 
over the crystal while holding the effect of the strong 
AP-coupling. The strong spin-coupling, however, real- 
izes only in the very narrow energy range of bandedge. 
The results for V(ui) and Q(u>), shown in Fig. ??, suggest 
that Q(uj) = —ui/A and V(u>) ~ in the wide range of 



energy oj as (—A < u> < A). This is consistent with the 
change in R(u>). As the energy w increases beyond —A, 
the probability of the carrier to reside Mn-site decreases 
so that the value of V(oj) rapidly approaches 0. 

In Fig. |3 the exchange-splitting energy, w p (up) — 
w p (down)[= —AE ex ], is plotted as a function of x(S z ). 
The data for each x are well fitted by a straight line. The 
slope corresponds to \N f3\ (or —AE ex /x(S z ) = —N f3). 
With the increase in x, the slope of the line decreases. It 
is worth noting that the straight line with x — 1 agrees 
with that of VCA. Thus, the present result suggests that 
the drastic decrease in \AE ex \/x{S z ) slope when x in- 
creases. 

Here, we compare the present result with that previ- 
ously obtained by the second order perturbation. Con- 
sidering the second-order scattering process due to the 
exchange interaction, Bhattacharjee 8] showed that not 
only the coefficient proportional to x{S z ) but also the co- 
efficient proportional to x(x(S z ))z is included in the term 

of (-^) . The latter coefficient comes from that the first 
order correction is taken in to account in the intermedi- 
ate process of the second-order scattering. According to 
his result, thus, the exchange splitting energy is not pro- 
portional to x(S z ). On the other hand, the present result 
suggests the proportional relationship between AE ex and 
x(S z ) holds in the wide range of the parameters IS, Em 
and x. The experimental observation in Zni-^Mn^Te 6| 
and Cdi-^Mn^Te shows that the spin-splitting AE as 
a function of x(S z ) is a straight line and that the slope 
No(a — 0) exhibits a large decrease for high x values, 
which seems to support the present result. 
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II-VI-based DMS 
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FIG. 4: Same as FigEJbut for IS/ A = -0.3 and E M /A 
0.0. 
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II-VI-based DMS 
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FIG. 5: Same as FigEJbut for IS/ A = -0.4 and E M /A 
+0.1. 



C. Some other results for II-VI-based DMSs 

Figures 01 and 0are the same results as Fig. Olbut for 
IS = -0.3A and E M = 0.0, and for IS = -0.4A and 
Em — +0.1 A, respectively. Figures indicate that the lin- 
ear relationship between exchange splitting and x(S z ) is 
kept for each case. Furthermore, the results show that 
not only the reduction in \IS\ but also positive Em sup- 
press the effect of multiple scattering. This can be un- 
derstood as follows. The repulsive interaction due to Em 
prevents the carrier to stay on Mn site longer. Thus, pos- 
itive Em substantially make the effect of the exchange 



interaction weaken. On the other hand, negative Em as- 
sists the carrier to reside on Mn site, resulting in appar- 
ently large |iVo/3|. Therefore, the present results request 
us reexamine the VCA and the hypothesis of Nq/3 = 21 
which have been widely accepted for II-VI-based DMSs. 
We need precise information on the exchange energy, the 
bandwidth and the band offset energy for more quanti- 
tative comparison between theory and experimental ob- 
servation. 



IV. RESULTS AND DISCUSSION FOR 
(GA,MN)AS 



A. Model parameters for Gai-^Mria; As 

In this section, we discuss in detail the reason why a 
positive Nq/3 is experimentally observed by magnetoop- 
tical measurement of ( Ga,Mn)As. For Gai^Mn^As, we 
take A = 2 eV [HHa, IS = -0.4A and E M = -0.3 A, 
as the same as previous work |l4j. These parameters 
lead to an impurity level at the energy of E a = — 1.057A 
in the dilute limit (x — ► 0), which is consistent with an 
acceptor energy of 0.113eV (= 0.057A) 28]. With the 
increase in x an impurity band forms around the ac- 
ceptor level. The impurity band merges into the host 
valence band at x > 0.035 when (S z ) — 0, while the 
down-spin DOSs unite at x > 0.017 when {S z ) = S. The 
result roughly agrees with the experimental observation 
of impurity-band- like states [2^, |3(J ■ Also this may be 
related to the insulator-metal transition 

reported to occur at x ~ 0.03 [3ll l32| . 



B. Case of low dilution 

As the typical case that x is so small that an impurity 
band forms separating from the host band irrespective of 
(S z ) = S, we investigate the case with x = 0.005; the 
results are shown in [H] ^ ??. As is shown in Fig. |6ja), 
the magnetic impurity band forms around the impurity 
level and imitates the model band. The numerical re- 
sult for the optical absorption spectrum A(oj) is shown 
in Fig. EJb). We notice that the peak of optical ab- 
sorption spectrum is almost near the bottom of the host 
band, although the impurity band exists in lower energy 
region. This can be explained as follows. The states of 
the impurity band are composed from the states of the 
wide range values in k space. The k = component in 
the impurity band states is nominal. Since the A(u>) is 
related to k = state, the optical absorption spectrum 
takes negligible values in the impurity band. 

As the consequence, the optical band edge, u> p , ex- 
ists almost at the bottom of the host band. Here, we 
should stress that the energy of the peak in A(u>) of 
up-spin carrier is lower than that of down-spin carrier 
(or Wp(up) < Wp(down)), although we have assumed the 
AFM exchange interaction (IS < 0). This implies that 
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the direction of the shift of the optical bandedge is op- 
posite from that predicted by the VCA. This may ex- 
plain the reason why the FM exchange interaction was 
reported on the basis of the polarized magnetoreflection 
measurement [9|, although the exchan ge i nteraction be- 
tween p holes and d spins is AFM |TolTlTl fl^| . 

In order to verify our picture, we first investigate the 
manner of coupling between the carrier spin and local- 
ized spins. We show the results for P(oj) and Q(w) in 
Fig.??. The result that V(u) » %^ and Q(u) w 1 in 
the impurity band indicates that the strong antiparal- 



lel spin coupling (AP-coupling) between carrier spin and 
localized spins occurs therein. Near the bottom of the 
host band, however, we find that V(oS) « SU S" 
gesting the parallel spin coupling (P-coupling) therein. 
On the other hand, Q{u>) takes high values suggesting 
the weak AP-coupling near the bottom of the host band. 
Only when as (S z ) > 0.955", Q{uj) takes negative values 
near the bottom of the host band. Here we should no- 
tice the difference in V(uj) and Q(ui). The optical carrier 
polarization, V{uo) 1 describes the spin coupling between 
the carrier's spin and localized spins averaged all over 
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the sites, while Q(uS) represents the strength of the spin- 
coupling at M-sites. Hence, the result suggests that the 
bottom of the host band apparently behaves as if the car- 
rier spin /erromagnetically couples to the localized spins 
although it may anii/erromagnetically couples to the lo- 
calized spins at M-sites. 

Next, we investigate how different way carrier states 
with up- and down-spin shift as the magnetization de- 
velops. The total number of states in the impurity band 
is x, irrespective of (5*2). When (S z ) = 0, the carrier 
states with up- and down-spin are completely the same, 
and the numbers of up- and down-spin states in the im- 
purity band are x/2, respectively. With increase in (S z ), 
as shown in Fig. EJa), the density of states with up-spin 
in the impurity band is suppressed, and finally vanishes 
when (S z ) = S. On the other hand, the number of states 
with down-spin in the impurity band increases from x/2 
to x. The total number of states keeps being 1.0 per a site 
for each spin (up and down) [see Eq. (|2.16|) ]. With the 
increase in (S z ), therefore, the carrier states with up-spin 
shift from the host band to the impurity band, whereas 
the carrier states with down-spin shift from the impurity 
band to the host band. 

To examine the mechanism in more detail, we show 
the result for local DOS in Fig. ??. Figure ??(c) indi- 
cates that rather large rate of the impurity band is com- 
posed from M-site states [R(u) w 0.5], in spite of small x 
(x = 0.005). The total number of impurity states, 0.005, 
is composed from M-site component of 0.00246 and A- 
site component of 0.00254, irrespective of (S z ). The re- 
sult shown in Fig. ??(a) indicates that the shift of the A- 
site states occur between the impurity band and near the 
bottom of the host band. Figure ??(b), on the other hand, 
indicates that the shift of the M-site states occur between 
the impurity band and over the wide energy range of the 
host band. This can be explained as follows. The effective 
local potential for carriers at M-site is Em + IS = — 0.7 A 
for AP-coupling while Em — IS = +0.1 A for P-coupling. 
The AP-coupling states constitute the magnetic impu- 
rity band, whereas the P-coupling states extend over 
wide range of host band due to the positive effective lo- 
cal potential. Thus, the up-spin carrier states at M-site 
are the AP-coupling sates in the impurity band when 
(S z ) — 0, whereas the P-coupling states spreading over 
the wide-range energy of host band when (S z ) — S. The 
carrier states with down-spin shift in the opposite way. 
Hence, nominal change occurs near the bottom of host 
band in M-site component DOS, as shown in ??(b). The 
present result reveals that the shift in the bottom of the 
host band is mainly ascribed to that in ^4-site compo- 
nent DOS. The spin dependent shift of the ^4-site states 
near the bottom of host band results in the shift of opti- 
cal bandedge. The direction of the shift is opposite from 
that predicted by the VCA. 

In Fig. we show the present result for the optical 
bandedge energies uj p with up- and down-spins as a func- 
tion of (S z )/S for x = 0.0005 and 0.005. The result sug- 
gests that the linear relationship between oj p and {S z )/S 



well holds in such dilute case that the impurity band 
forms separate from the host band. In Fig. [8J we plot 
the exchange splitting cj p (down)— oj p (up) [= +AE ex ] as a 
function of x(S z )/S for x = 0.0005, 0.001, and x = 0.005. 
The data are well fitted by a straight line. The Nq/3 
with S — 5/2 deduced from the slop of the straight 
line is +1.31 eV. Szczytko et. al measured the exciton 
splitting in Gai-^Mn^As with x = 0.00047, 0.00027 and 
0.00022 by polarized magnetoreflection, and showed that 
the data are proportional to magnetization, to obtain 
Ao/3 = +2.5 ± 0.8 eV @. The agreement between the 
present result and experimentally obtained one is satis- 
factory. Therefore, we conclude that the spin-dependent 
shift of the carrier states between the impurity band and 
host band accompanying with the change of magnetiza- 
tion causes the apparently FM behavior of the optical 
bandedge. 
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C. Case of moderate dilution 

As the typical case that x is so large that the impurity 
band merges to the host band irrespective of {S z }, we 
investigate the case with x = 0.04; the results are shown 
in FigEI ~ ??. In contrast to the case of x = 0.005, 
the impurity band and host band have already united to 
form the band tail, as shown in Fig[!Ja). The peak of 
the optical absorption spectrum is broaden due to the 
increase of the disorder which is caused by the additional 
corporation of M ion, as shown in Fig^b). However, the 



peak of optical absorption spectrum is not in the bandtail 
originated from the impurity band but still exists in the 
energy range near the bottom of the original host band. 
Thus, the direction of the shift in the optical band edge 
is opposite to that predicted by the VCA, as is the case 
of x — 0.005. The results for V(ui) and Q(u>), shown in 
Fig.??, indicate that the strong AP-coupling realizes in 
the bandtail, whereas Viw) shows the weak P-coupling 
at the energies near the optical bandedge. The result for 
the shift of the states accompanying the change in (S z ), 
deduced from Fig ??, has similar tendency with that of 
x = 0.0005, although the mergence of the impurity band 
makes the feature less clear. 

In Fig llOl we show the optical bandedge energies, 
Wp(up) and w p (down), as a function of (S z )/S for x = 
0.02 and x = 0.04. When x = 0.04, the absorption 
spectrum of up-spin is broadened. The exchange split- 
ting AE ex increases monotonically with the increase in 
(S z ) I S, while the linear relationship hardly holds. In Fig. 
ITDwe plot the exchange splitting, w p (down) — w p (up) [= 
+AE ex ], as a function oix{S z )/S for x = 0.02,0.03, and 
0.04. The data do not fit a straight line well. Note that 
for 0.017 < x < 0.035, the impurity band separating 
from host band when (S z ) — unites to the host band 
when (S z ) = S. The Nq/3 deduced from the slope of the 
straight line in Fig. ^]is +1.36 eV, which is consistent 
with the experimental observation (see below). 

Here we compare the present result with that of mag- 
netoabsorption fl2|. Note that the band-gap energy of 
GaAs is 1.52 eV [Tl|. In the absorption spectrum (see 
Fig. 1 in Ref. 12), a rather weak structure is visible be- 
low 1.5 eV. We regard the structure as the optical transi- 
tion related to the impurity band, although Szczytko et 
al. disregarded it as the below-the-gap transition. When 
x is so large that the impurity band units the host band, 
the optical absorption spectrum A(u>) takes appreciable 
values at the energies in the impurity band, as shown 
ininib). In the energy range above 1.5 eV, on the con- 
trary, the absorption band increases monotonously with 
the increase in the photon energy, which we assign to the 
optical absorption related with the host band. Under 
the external magnetic field, the edge splits about 0.1 eV, 
which is the feature of sp-d exchange effect but opposite 
way than for Cdi-^Mn^Te. Evaluating the spin splitting 
as the relative edge shift at high energies as 1.7 eV, they 
obtained N {a - 0) = -2.1 eV for x = 0.032 and -1.7 
eV for x = 0.042. Thus, they measured the energy shift 
at the host band to obtain the No(a — (3). The result is 
consistent with our view that the optical bandedge near 
the bottom of original host band behaves as if the ex- 
change interaction is FM although the AFM exchange 
interaction acutely operates between the carrier and lo- 
calized spins at Mn site. The values of No(a — (3) are 
comparable with No(a — (3) = —(2.3 ± 0.6) eV obtained 
in the very dilute cases of x = 0.00022 - 0.00047 0. 
Szczytko et al. concluded that the coincident of values 
of No(a — (3) for different x is accidental because they 
believed the splitting inversion is so caused by the Moss- 
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Burstein effect that the splitting depends on the carrier 
concentration. The coincidence, however, might not be 
accidental because the spin splitting near the bottom of 
the host band was measured in both cases. The states 
near the bottom of the (original) host band are always al- 
most empty, although the carriers may enter into the im- 
purity band or the bandtail originating from the impurity 
band. Since the optical transition is related to the states 
near the bottom of the (original) host band, as shown 
in this study, the carrier concentration would not signifi- 
cantly affect the exchange splitting optically observed in 
III-V DMSs. The impurity band and/or bandtail struc- 
ture changes with the Mn concentration x. In the ab- 
sorption band of Gai-^Mn^As experimentally obtained, 
the impurity-like structure of x — 0.042 more overlaps 
with the host-band-like structure than that of x — 0.032. 
This may be another supporter of our picture. 



the impurity band forms separate from host band, the 
optical bandedge exists not at the bandedge of the 
impurity band but near the bandedge of host band, as 
shown in this study. The optical bandedge behaves as if 
the exchange interaction is ferromagnetic although the 
antiferromagnetic exchange interaction actually operates 
at Mn site. We concluded that the spin-dependent shift 
of the carrier states between the impurity band and 
host band causes the apparently ferromagnetic shift 
of the optical bandedge. The conclusion is valid even 
when x is so large that the impurity band unites the 
host band because the optical bandedge exists near the 
bottom of the original mother band. The present result 
is consistent with the theory for the mechanism of the 
carrier-induced ferromagnetism in III-V based DMS's 
that we have previously proposed |14j . 



V. CONCLUDING REMARKS 

In this study, applying the dynamical CPA to a simple 
model, we investigated the behavior of the optical band- 
edge in DMSs in a systematic way. For AY_ x Mn x B Yl - 
type DMS, the present study reveals that the linear re- 
lationship between exchange-spitting AE ex and the av- 
eraged magnetization ^(S^)! widely holds for different 
values of x. The ratio, AE ex /x(S z ), however, depends 
not only the exchange energy (IS) but also the band off- 
set (Em)- The present theory can qualitatively explain 
the x dependency of Nnta — j3) reported in Zni-^Mn^Te 
@ and Cdi_ x Mn x Te 0. For the quantitative descrip- 
tion, however, we need more precise knowledge on the 
exchange energy, the bandwidth, and the band offset en- 
ergy. This will be discussed in detail elsewhere. Here, 
we have to indicate that the present model does not 
take into account of many features such as multiband 
effects, band anisotropy, and excitonic structure, which 
exist in real DMS's. In the CPA, furthermore, the col- 
lective mode, correlation and/or clustering effect of lo- 
calized spins, which may become significant for large x 
region, are completely out of scope. These issues remain 
for future study. 

Regarding the p-d exchange interaction of 
Gai-xMn^As, there have been long controversial 
discussions not only on the amplitude but also even on 
the sign. The uncertainty has made the model for the 
carrier-induced ferromagnetism difficult to establish. In 
the present paper, we have proposed a new interpre- 
tation for the experimental result in magnetooptical 
measurements. A Mn 2+ ion in GaAs acts as both 
an acceptor and a magnetic impurity. Therefore, on 
the Mn-site, a carrier (p hole) is subject to the local 
potential which includes the p-d exchange interaction 
together with the attractive Coulomb potential. In 
the dilution limit, thus, there appears an acceptor 
level. With the increase in x, an impurity band forms 
around the impurity level. In such low dilution that 
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APPENDIX A: DYNAMICAL CPA — t MATRIX 
FORMALISM 

1. t matrix elements of A ion embedded in the 
effective medium 

Here we omit the site suffix. The t matrix elements, 
which represent the multiple scattering of carriers with 
| (j)-spin due to the A ion potential Ea embedded in 
the effective medium Ef (Ej,), is given by (l(i| 



t A - 



1-(£U-£ T )F T ' 
1-(E A - Ex)^ • 



(Ala) 
(Alb) 



Here, F^cv)] is the diagonal matrix element of a 

propagator P with respect to the effective medium, E^ [= 
E M (w)] (fi =| or |), and is calculated by 



F M = ( M |P| M ) = f de 

J-A 



P(e) 



LU - £ - E^U;) 



, (A2) 



where p(to) is the model DOS. For p(tu) given by Eq. 
I|2.14|) . we obtain 



F»A = 2 >< U ^ 



Note that tf t = = 0. 



'u - E, 



- 1 XA3) 
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2. t matrix elements of M ion embedded in the 
effective medium 



Then, Eq. G3 x (F^)"* + Eq. C5| x^-F^;)- 1 
leads to an equation including tff only (^ is canceled): 



In accordance with the definition of the t matrix, Eq. 
(|2.7|) . we have 



t M [l-Pv M ] = v M . 



(A4) 



Equation (|A4(I is written in the spin-matrix element ex- 
pression as 



t 



M + M TP ,, M + M P *, M 
TT _ Mt^T^fr ~ hi*l v ll 



TT > 



TT — Ti J- U — TT T w Ti ~ TT 



(A5) 
(A6) 



(A7) 



By using the following definitions and/or symbols intro- 
duced for simplicity, 



by a similar procedure. The resulting expressions are 



Vi 

v M 

TT 

,,M 

Ui 



; T1 



T > 



,M 



E M - IS Z 
vn = E M + IS Z - S| , 
-IS- 
-IS+ 

E M - I(S 2 - 1) - E T , 
F M + + 1) - Ej, , 
7 2 S'_5 + = I 2 [S(S+1) 
I 2 S+S- = I 2 [S(S + 1) 



Si 
Si 



s z 
s z 



(A8a) 
(A8b) 
(A8c) 

(A8d) 
(A8e) 
(A8f) 
(A8g) 
(A8h) 



and recalling the commutation relationships between the 
components of S, 



SS Z = (S z + 1)S- , 

(s+r 1 = [s(s + i)-(s z f 



(A9) 

S.^OSL) ,(A10) 



V r + F l (W 1 -V 1 U i ) 



1-F^-F^+F^-V^)] 

vx+n^T-n^T) 

1 - FfC/f - F;[Vj + F T (^ X - y^)] 
1 

1 - F l U i - F T [V T + F X (W T - VtC/x)] 
(-75_) ' 



(Alia) 
(Allb) 



(-LS_) 



,M _ 

TT - 



1 - FffTf - F ; [^ + F T (Wj - y x (7 T )] ' 

(Allc) 

1 



1 - F T £/ T - F ; [^ + F T (W X - V1C/ T )] 
= {-IS. ' ' 



(-LS 4 



1 - F l U l - F T [V T + F^W, - V^)] ' 

(Alld) 

V| (V|) is the spin-diagonal component of the interaction 
between a carrier with j (J.) spin and the local potential 
on the M ion embedded in the medium of Ej (EjJ. A 
carrier with f (I) spin which has already flipped in the 
previous scattering is subjected to C/| (f/jj on the M ion 
embedded in the medium, wherein the d spin operator 
S z is replaced by S x — 1 (5 Z + 1). Furthermore, Wf (Wj,) 
represents the interaction energy required by a carrier 
with | (!) spin to flip and then reverse its spin. 



APPENDIX B: DYNAMICAL CPA 
FORMALISM 



LOCATOR 



CPA locator condition 



we obtain an explicit expression for using no more 
approximations. Other t- matrix elements are obtained 



In this subsection, we briefly outline an alternative but 
equivalent condition of the CPA, which is an extension 
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of the CPA using the locator formalism [33j. Assum- 
ing that the spin-dependent effective medium surrounds 
an arbitrary site n, we consider the transfer of the car- 
rier with spin \i between site n and the effective medium 
(i.e., S| and £jj by the site-renormalizcd intcractor J M . 
Then, the propagator G A (G M ) associated with the real 
potential of u A (it„ ) embedded at site n in the medium 
is defined by 



G' 



1 



L n S fj, J r ti a nfiO>nn 



G 



M 



(Bl) 
(B2) 



When we set the coherent potential £ M on the site n in 
the effective medium, the reference propagator, 



P 



1 



(B3) 



is equivalent to the Green function for the effective 
medium. Thus, the diagonal matrix element of P is equal 



to defined by Eq. (HJ|: 



F^w) = (rt[i\P\nfi) = 



(B4) 



Equation i|B4() gives the relationship between and F M ; 
£ M = 1/(oj — S M ) is called a locator. Hereafter, for 
the sake of simplicity, the site-diagonal elements in the 
Wannier representation (nfi\G \nv) are written as G A V 
(fJ,,v =T, or |). Then, the spin-diagonal element of G A 
is given by 



F A (ui) 



1 



F, 



Ju 



(B5) 



and the spin-off-diagonal elements are G A ^ = G A ^ = 0. 
The site-diagonal elements of G M are obtained after a 
somewhat complicated calculation using the commuta- 
tion relationships between the components of S but with 
no further approximations, as 
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+ S Z ] 


[uj - (E M 
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[uj - (E M + IS Z ) - J t ][uj - E M + 7(5, - 1) - J T ] - P[S(S + 1) - 5| + S z ] 



(-IS-) 



(-75+) 



(B6a) 
(B6b) 

(B6c) 
(B6d) 



Note that the site-diagonal elements of G M involve spin 
operators. Thus, F^ 1 (uj) is defined as the thermal aver- 
age of the spin-diagonal element Gff^ by 

S s 
f» = C= E Gf M (5,)exp(A5,)/ £ exp(A5 z 

s z =-s s z =-s 

(B7) 

where A (= h/ksT) is determined so as to reproduce a 
given value of (S z ) [see Eq. (|2.13|) ]. Note that the spin- 
off-diagonal elements (Gf[) = (G$) = 0, because Gf{ 
(GH) includes 5_ (5+) in their final form. Finally, the 
CPA condition in the locator formula is given by 

7» = (1-x)f£(w)+xFF(u). (B8) 

When F^ is given, J M is calculated by Eqs. I|A2() and 
(|B4)l [i.e., = (A 2 /4)F M for the model band defined by 



Eq. (EH]. Then, F A and F]f are calculated by Eqs. 
(|B5|) and (|B7|) . and consequently, F^ is again obtained 
by Eq. ljB8|) . Therefore, F M and J M are determined self- 
consistently. 

2. Local DOS 

The advantage of the locator formula CPA is that it is 
straightforward to determine species-resolved DOS, i.e., 
DOS associated with each kind of ion in the alloy. The 
local density of states (local DOS) at the A (M) site is 
obtained by 

7>» = --In<M , (B9a) 
Dff(u) = -ilmTfH . (B9b) 
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In actual calculations, first we determined and £ M 
by the t matrix formula CPA, and then calculated J M , 
and consequently D^(u>) and Djf(cu), We numerically 
verified the relations 

/OO POO 
Dfa)du= / Dff(u)du = l, (BIO) 
-oo J — OO 

and 

D^) = (l-x)Dfa)+xD?f(u) , (Bll) 



which are a consequence of the CPA locator condition, 
Eq. CEBJl. 



3. Spin-coupling strength Q(cu) 

We define the spin-coupling strength Q(u>) as the nor- 
malized (negative) inner product between the carrier spin 
and localized spin on the M site: where D M (to) = 
Df{u) + Df(uj) and 



CH = ^2(S(u - H) a^a^u ■ Saw) 

- llm'^(nv\G h \^)\n^)(n^\al lti cr^ ■ Sa nv \nv) 



--Im((G^ -G A L l)S z + Gf{S+ + G A L \S-) . (B12) 
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